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The  theory  of  reliability  estimation  developed  in 
previous  reports  (WADD  TR61-53,  ML-TDR-64-300)  has  been 
applied  to  three  types  of  aircraft,  a  civilian  transport, 
a  heavy  bomber  and  a  fighter  designed  by  current  proce¬ 
dures,  for  \tiiich  operational  records,  multiple  structural 
tests  and  records  of  service  experience  are  available. 
Failure  rates  for  critical  ultimate  load  conditions  have 
been  evaluated  on  the  basis  of  data  obtained  from  various 
sources  and  compared  with  service  experience.  Lives  as¬ 
sociated  with  equal  risk  of  ultimate  load  failure  and  fa¬ 
tigue  failure  (or  initial  structural  fatigue  damage)  have 
also  been  computed. 

The  obtained  numerical  values  vrtiich  reflect  cur¬ 
rent  design  practices  can  serve  as  the  basis  for  a  ra¬ 
tional  comparative  reliability  analysis  of  new  designs 
involving  new  materials  and  different  design  criteria 
and  missions  spectra  and  profiles. 
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The  theory  developed  previously  for  determining  the 
probability  of  structural  failure  under  both  ultimate  and 
fatigue  loading^  has  been  applied  to  certain  aircraft  struc¬ 
tures  for  which  operational  records  and  multiple  structural 
tests  are  available  and  the  results  of  the  analysis  have 
been  compared  with  results  from  service  experience.  Jin 
analysis  has  also  been  made  of  the  effect  of  various  design 
and  specification  parameters  on  the  probability  of  failure 
under  ultimate  load. 

The  method  used  takes  account  of  the  statistical 
variation  of  operational  loads  acting  on  the  aircraft  as 
well  as  of  the  ultimate  strength  and  fatigue  life  of  the 
structure  and  has  been  applied  to  the  following  representa¬ 
tive  cases:  (a)  civil  transport  aircraft,  (b)  heavy  bomber 
aircraft,  and  (c)  fighter  aircraft.  In  each  case  the  "risk" 

of  ultimate  load  failure  (r  ) ,  the  risk  of  fatigue  fail- 

1  rF 

ure  (r  )  and  the  fatigue  sensitivity  factor  [f(N)  *  —  ] 
r  r 

u 

have  been  determined. 

For  ultimate  load  failure  reasonable  agreement  with 
service  experience  has  been  found  and  the  maximum  values  of 
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the  confuted  failure  rates  associated  with  current  design 
procedures  are  as  follows: 

..7 

(a)  Civil  transport:  10  per  hour  (thunderstorm 

downgust  critical) 

—8 

(b)  Heavy  bomber:  10  per  hour  (lateral  gust  on 

tail  assembly  critical) 

-5 

(c)  Fighter:  10  per  hour  (upward  maneuver 

load  critical) 

The  actual  failure  rates  may  vary  by  several  orders  of  magni¬ 
tude  depending  on  the  assumed  design  and  operating  conditions 
and  the  proposed  method  of  reliability  analysis  makes  it  pos¬ 
sible  to  evaluate  the  effect  of  a  change  in  such  conditions 
on  failure  rate .  Comparison  with  service  experience  indi¬ 
cates  that  the  theory  provides  an  adequate  method  of  relia¬ 
bility  prediction. 

The  risk  of  fatigue  failure  is  a  function  of  the 
life  (N) ;  introducing  the  life  Nq  at  which  the  risks  of 
ultimate  load  and  fatigue  failure  axe  equal  as  a  design  cri¬ 
terion  these  lives  Nq  at  which  for  current  design  proce¬ 
dures  the  structures  become  fatigue  critical  are : 

(a)  Civil  transport  (design  service  life  for  ultimate 

load  assumed  to  be  20,000  hours) 

N  *  15,000  hours 
o 

(b)  Heavy  bomber  for  mission  profiles  including  low 
level  flight  (initial  fatigue  failure) 

N  *  960  hours 
o 

(c)  Fighter 

N  *=  1,900  hours 
o 
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The  associated  reliability  figures  for  ultimate  load  design 
are  of  the  order  of  R  *  0.99  for  (a)  and  (b)  and  R  *  0.94 
for  (c) .  The  above  figures  indicate  that  for  currently  as¬ 
sumed  operational  lives  and  considered  operational  load 
spectva  the  risk  of  fatigue  failure  is  critical  in  all  cases 
considered  in  this  analysis. 

Having  established  the  numerical  values  for  ultimate 
load  failure  rates  and  for  fatigue  lives  at  equal  risk  of 
ultimate  load  and  fatigue  failure  reflecting  current  design 
practices,  a  rational  comparative  reliability  analysis  of 
new  designs  becomes  possible  using  different  design  cri¬ 
teria,  different  structural  materials  and  different  mission 
spectra  and  profiles,  but  retaining  the  method  of  analysis 
outlined  and  illustrated  in  this  report. 
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1.  INTRODUCTION 


The  rising  importance,  in  recent  years,  of  the  prob¬ 
lem  of  structural  reliability  is  largely  a  result  of  the  in¬ 
creasing  size  and  complexity  of  aircraft  structures  as  well 
as  the  rapid  growth  of  aircraft  operations . 

To  insure  adequate  safety  of  civil  transport  air¬ 
craft,  it  is  necessary  that  the  probability  of  structural 
failure  be  negligible  compared  to  the  normal  hazards  in¬ 
volved,  as  in  any  form  of  transport.  In  the  case  of  mili¬ 
tary  aircraft  the  risk  of  structural  collapse  should  be  much 
lower  than  the  accepted  operational  risk.  To  meet  these 
requirements  the  probability  of  structural  failure  at  any 

stage  must  be  determined.  While  this  problem  has  attracted 

2  3  4 

the  attention  of  a  number  of  authorities,  *  a  solution 
of  the  general  case  in  terms  of  known  parameters  has  not 
yec  evolved. 

An  approach  to  the  problem  has  recently  been  made 
based  on  reliability  theory,^-  and  the  purpose  of  the  present 
investigation  is  as  follows: 

(a)  to  investigate  the  probability  of  fatigue  and  ul¬ 
timate  load  failure  when  the  various  design  and  operational 
parameters  are  varied; 
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(b)  to  use  data  obtained  from  service  operations 
and  to  compare  the  performance  predicted  with  that  actually 
achieved  for  both  civil  and  military  aircraft  types. 

These  objectives  have  required  the  assembling  of  com¬ 
prehensive  data  on  the  ultimate  strength  and  fatigue  strength 
of  various  structures,  as  well  as  flight  load  data  from  a 
variety  of  aircraft  operations  and  service  data  from  civil 
and  military  aircraft  including  records  of  structural  fail¬ 
ures  . 

The  program  has  been  supported  by  the  Research  and 
Technical  Division  of  the  United  States  Air  Force  Systems 
Command ;  the  continuous  interest  and  active  assistance  in 
all  phases  of  this  program  of  Mr.  W.  J.  Trapp,  Chief,  Strength 
and  Dynamics  Branch,  Air  Force  Materials  Laboratory,  is  grate¬ 
fully  acknowledged.  The  program  has  received  cooperation 
from  the  Aeronautical  Research  Laboratories  in  Australia; 
the  Research  Institute,  University  of  Dayton;  the  United 
States  Naval  Air  Engineering  Center;  the  National  Aeronau-  \ 

tics  and  Space  Administration;  the  Federal  Aviation  Agency; 
leading  aircraft  manufacturing  firms  in  the  United  States; 
the  Royal  Aircraft  Establishment  in  the  United  Kingdom;  and 

i 

1 

Sud-Aviation  Aircraft  Company  in  France.  i 

i 

i 

i 
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2.  FLIGHT  LOADS 

Flight  load  data  have  been  obtained  from  various 
sources  and  have  been  presented  here  in  two  braod  divisions: 
atmospheric  gusts  and  aircraft  maneuvers.  All  data 
available  were  in  the  form  of  frequency  distributions  which, 
for  our  purposes,  have  been  transformed  into  probability  dis¬ 
tributions  associated  with  a  value  for  the  total  number  of 
load  occurrences  per  hour  (or  per  mile).  Equations  for  the 
probability  distributions  and  the  total  frequencies  of  oc¬ 
currence  appear  in  Table  1. 

(a)  Atmospheric  Gusts 

These  data  refer  primarily  to  the  effect  of  atmos¬ 
pheric  turbulence  ;  the  effect  of  any  maneuver  loads  that  occur  is 
included  since  these  effects  are  minor  where  gust  loads  are  relevant. 

(i)  Thunderstorm  Gust  Data 

Results  from  the  National  Advisory  Committee  for 
Aeronautics  investigations  in  1941-42,  1946  and  1947  on  gust 
frequencies  in  thunderstorms  have  been  obtained  from  Tolef- 
son.^  A  regression  analysis  has  been  carried  out  and  the 
various  results  have  been  combined  to  give  the  probability 


spectrum  presented  in  Fig.  1  associated  with  an  estimated 
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14  gusts  per  mile.  The  average  percentage  (a  %)  of  the 

total  distance  spent  in  thunderstorms  has  been  estimated 

as  a  function  of  operating  height these  data  have  been 

reproduced  in  Table  2.  For  an  average  cruising  speed  V 

c 

the  average  number  of  loads  per  hour  is  therefore 
14  x  2  x  10“2  x  Vc  *  0.14a  Vc. 

(ii)  General  Datafor  Civil  Transports 

Comprehensive  gust  load  data  have  been  obtained 

by  the  Royal  Aircraft  Establishment  from  civil  transports 

7 

operating  over  world-wide  routes.  These  data  are  repro¬ 
duced  as  a  probability  distribution  in  Fig.  2.  The  total 
nunber  of  gusts  per  nautical  mile  is  shown  as  a  function  of 
height  in  Fig .  3 . 

(iii)  Lateral  Gust  Data 

The  only  available  data  on  lateral  gust  frequen¬ 
cies  aure  those  recently  obtained  by  the  University  of  Dayton 

8 

on  United  States  heavy  bombers .  Data  for  both  high-  and 
low-level  operations  have  been  presented  as  probability  dis¬ 
tributions  in  Fig .  4  with  the  associated  gust  counts  per  hour . 


(b)  Maneuver  Accelerations 
(i)  Fighter -Aircraft 

Maneuver  load  data  on  fighter  aircraft  from  a 
variety  of  sources  are  reproduced  in  Fig.  5.  The  United 
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States  and  British  data  are  the  average  of  a  large  number  of 
records  on  different  fighter  aircraft  types .  The  United 
States  data  are  based  on  fighters  in  the  Air  Force  FII  cate¬ 
gory  with  an  ultimate  load  factor  n  =11. 

u 

Variation  occurs  according  to  role,  the  ground-attack 
mission  giving  more  frequent  positive  (upward)  acceleration 
counts  than  any  other  mission.  However,  the  distribution 
with  the  density 

.  2 
_  An 

2  2a  3 

p(An)dAn  =  /Tf  °n  6  dAn  An  >  °  (2*1) 

has  been  found  to  give  a  consistently  good  representation  for 
the  positive  accelerations  An  except  at  very  small  accelera¬ 
tion  values.  Eq.  (2.1)  represents  the  density  function  of  a 
normal  distribution  with  mean  zero  truncated  at  the  mean  for 
positive  values  with  normalizing  factor  2.  All  missions  can 
therefore  be  represented  by  taking  suitable  values  for  the  mission 

parameter  a  as  shown  in  Table  I . 
n 

This  procedure  is  quite  satisfactory  for  the  predic¬ 
tion  of  ultimate  load  failure.  For  general  application  in 
fatigue  calculations,  however,  it  is  recommended  that  an 
exponential  distribution  be  added  for  both  positive  and  nega¬ 
tive  loads  associated  with  their  total  occurrences  per  hour 
to  give  a  frequency  spectrum  of  fatigue  loads . 
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(ii)  Bomber  Aircraft 


Combined  gust  and  maneuver  load  data  on  United 
States  heavy  bombers  are  presented  as  a  probability  distri¬ 
bution  of  accelerations  in  Fig.  6.  Unlike  the  preceding 
data,  which  are  of  general  application,  these  results  refer 
to  the  particular  type  only.  In  order  to  facilitate  their 
application  in  the  calculations  in  Section  3,  these  data  are 
presented  in  a  form  which  includes  the  effect  of  the  fleet 
operating  characteristics. 


3.  DISTRIBUTION  OF  STATIC  STRENGTH 


To  investigate  the  variability  of  ultimate  strength 
of  aircraft  structures,  nultiple-test  data  from  11  different 
types  of  structure  and  9  types  of  mainplane  panels  were  assem¬ 
bled  for  analysis  as  shown  in  Table  3.  The  panels  were  all 
loaded  to  failure  in  compression  but  the  different  types  of 
structures  were  tested  under  various  loading  cases  and  ten¬ 
sion,  compression  and  shear  failures  are  represented  in  these 
results.  The  data  from  these  19  groups  therefore  give  a  good 
representation  of  the  ultimate  strength  behavior  of  air¬ 
craft  structures  in  general . 

(a)  Analysis  of  Data 

Analysis  of  the  data  suggests  that  a  single  distri¬ 
bution  of  the  ultimate  failing  load  about  the  mean  value  can 
be  obtained,  irrespective  of  the  type  of  structure  tested  or 
its  mode  of  failure. 

To  test  this  hypothesis  against  the  experimental  data, 
the  population  means  \±  have  to  be  estimated  in  each  group  by 

•K  • 

i 

the  experimental  mean  fL  which  for  some  groups  is  determined 
on  only  3  or  4  specimens.  In  an  attempt  to  reduce  this  bias,  the 
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variate  "£  =  —  has  been  represented  by  the  standardized 
ax 

var iate 


R.  -  R.  R.  -  R. 

11  11 


which  has  been  estimated  by 


Y 
m 

The  values  R^SX  for  estimating  or  are  based  on  the  sample  means  R^ 
and  the  combined  standard  deviation  S  determined  on  the 

A. 

total  of  170  values, 
vraluesj 

Thej^T^ have  been  calculated  for  every  specimen  in 
each  group  and  then  arranged  in  ascending  order  of  magnitude 
(m  =  l,2,...,n).  These  values  have  been  plotted  against 
the  mean  relative  frequency  in  Fig.  7;  it  is  found 

that  the  n  =  170  points  so  obtained  are  all  distributed 
along  a  smooth  curve  with  no  tendency  to  segregate  according 
to  groups.  To  show  this  the  10  most  extreme  points  in  both 
tails  of  the  distribution  have  been  identified  by  their  group 
letter  (as  listed  in  Table  3)  and  the  extreme  points  for  each 
group  are  also  identified  in  the  array  of  170  points  by  their 
group  letter  . 
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The  probability  distribution  as  represented  by  the 
data  points  is  symmetric;  it  follows  the  normal  distribution 
closely  near  the  mean  but  diverges  from  it  considerably  in 
the  vicinity  of  the  extreme  values. 

Of  the  many  distributions  tried,  the  "t"  distribution 
was  the  only  one  giving  a  satisfactory  representation;  the  ex¬ 
cellent  fit  shown  in  Fig.  7  was  obtained  using  the  "t "  dis- 

3 

tribution  with  3  degrees  of  freedom  and  setting  Y  =  ~  t. 

Hence 


r . 


R  — R 

Vy)  =pyC  ■ 


i  x 


P  f 

P(vt=3)  V3  V 


(3.1) 


where  P  (t)  is  the  probability  function 
Vt  S 

tribution  with  3  degrees  of  freedom. 

X 

It  then  follows  that  the  variate 

SX 


of  the  "t"  dis- 


Ri~|iRi 

^R.SX 


based 


on  the  true  means 


has  the  distribution 


vt=3  vv3 


JL> 


(3.2) 


j 


! 
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while  the  true  standardized  variable  f  *  — —  must  have 

Vx 

the  distribution 


P(?)  «  P(/3  ?)  =  P  A/3  ~ 

Vt=3  ax 


(3.3) 


since  P  ,(t)  has  variance  of  3. 
vt=3 


in 


Equivalence  of  (3.3)  and  (3.2)  gives  a  -  ’s*rst’  S 

X  4  X 


1.299  Sx  for  best  fit  of  the  "t"  distribution  shown  in 


Fig .  7 . 

2  2 

A  comparison  of  both  S  and  a  with  the  sample 

A  A 

variances  of  the  various  groups  of  data  has  been  carried  out 

using  the  F  test.  The  results  presented  in  Table  3  show  that 
2  2 

both  Sx  and  ox  as  population  variances  give  no  signifi¬ 
cant  difference  from  the  sanple  variances  in  10  cases  out  of 
the  19.  However,  since  the  F  test  refers  to  the  normal  dis¬ 
tribution  it  will  prove  a  too  severe  test  for  the  compati¬ 
bility  of  variances  in  a  "t"  distribution  of  three  degrees 
of  freedom.  The  9  cases  in  which  there  is  a  significant 
difference  do  not  therefore  necessarily  reject  the  hypothe¬ 
sis  of  a  common  "t"  distribution.  But  the  results  do  indi¬ 
cate  that  comparison  with  the  sample  variances  does  not 

2 

reject  the  value  of  a  as  the  population  variance  in 

A 

2 

favor  of  the  value  of  S  . 

A 
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2 


is  therefore  taken  for  the  variance 


The  value  a. 


of  x  = 


R-m-t 


.  Hence 


P  (x)  =  P  U 

A 


of  =  V1+^)'2  = 

ox  tt  v  >  ctv  ,  — 

X  X  TUC  *  + 


X2) 


(3.4) 


Px(X)  =  f 

—  CJD 


,  .  .  V  .  1  X  .  1 

p(x)ctx  -  - - r“7 — rr  +  —  artan  —  +  -r 

n(ax2  +  x2)  tt  aX  2 


(3.5) 


Within  the  region  covered  by  the  data  (-.15  £  x  £  .15)  it 


has  therefore  been  possible  to  establish  for  the  variation 


in  relative  strength  x  = 


Mr 


a  unique  distribution  with 


a  known  variance . 


(b)  Generalized  Probability  Distribution 

The  validity  of  extrapolating  the  "t"  distribution 
beyond  x  =  0.15  is  of  little  practical  importance  since 
under  any  of  the  service  load  spectra  the  probability  of 
failure  in  the  region  x  >  0.15  is  negligible. 

For  the  critical  range  -l.C  £  x  <  -0.15,  however, 
extrapolation  of  the  "t"  distribution  proves  unacceptable. 
The  probability  density  function  is  finite  for  structures 
of  zero  strength  (x  =  -1.0),  which  is  physically  unrealis¬ 
tic;  it  predicts,  moreover,  a  high  failure  rate  for  struc¬ 
tures  of  negligible  ultimate  strength  when  combined  with 
typical  load  spectra . 
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This  is  referred  to  farther  in  Section  4  where  it 


is  shown  that  with  the  exponential  gust  load  spectrum  a  "t" 
distribution  of  structural  resistance  produces  a  spurious 
increase  in  the  marginal  density  function  of  structural 
resistance 

P(R)  =  P„(R)  [1  -  PS(R) J 

for  values  of  structural  resistance  less  than  80%  of  the 

,  S 

mean  va lue  (  x  =  - <  -0.2  ). 

This  problem  can  be  overcome  by  truncating  the  Mt" 
distribution  of  structural  resistance  at  a  suitable  value, 
on  the  quite  justifiable  assumption  that  structures  of 
lower  strength  will  be  eliminated  by  inspection  during  manu¬ 
facture  .  It  is  known  from  experience  that  the  rare  occur¬ 
rences  of  ultimate  failure  in  service  do  not  arise  because 
of  structures  having  extraordinarily  low  strength  but  rather 
because  structures  of  adequate  strength  encounter  extreme 
loading  conditions. 

Fundamentally,  however,  it  is  desirable  to  apply  the 
reliability  approach  over  the  whole  range  of  possible  struc¬ 
tural  resistance  and  this  also  makes  it  possible  to  assess 
the  various  parameters. 
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However,  until  further  data  on  ultimate  strength 

are  available,  the  behavior  of  the  distribution  near  the 

extreme  values  can  be  based  only  on  physical  reasoning . 

The  probability  distribution  of  the  ultinate  stress 

at  failure  in  metals  can  be  approximated  by  a  log-normal 
9 

distribution  or,  in  fact,  for  the  low  variances  that  apply, 
a  simple  normal  distribution.  The  ultimate  failure  of  struc¬ 
tures  is  further  dependent  on  the  load  eccentricity  in  the 

2  2  2  2 

critical  members,  represented  by  e  =  z^  +  z2  +  z^ 

where  z^,z2,z^  are  the  components  of  eccentricity  parallel 

to  any  three  fixed  perpendicular  axes.  If  the  deviations 

zl,z2,z3  are  assumed  to  be  normally  distributed  about  zero, 

2  2 

e  is  therefore  an  %  variate  with  three  degrees  of  freedom. 

This  gives  some  physical  explanation  of  why  the  struc¬ 
tural  resistance  appears  as  a  "t"  statistic,  defined  by  a 

normal  variate  (the  distribution  of  ultimate  stress)  divided 

2 

by  the  square  root  of  an  X.  variate  with  3  degrees  of  free¬ 
dom  (the  distribution  of  load  eccentricity  "e").  It  repre¬ 
sents  the  effect  of  load  eccentricity,  inducing  bending 
stresses  which  add  to  (or  subtract  from)  the  direct  stress 
in  members  at  the  critical  point. 

The  allowable  eccentricities  due  to  lack  of  straight¬ 
ness  in  members,  inexact  location  of  bolt  and  rivet  holes. 


-16- 


etc.,  will  be  held  within  fixed  limits  by  the  normal  produc¬ 
tion  standard.  Beyond  these  limits  it  seems  reasonable  to 
assume  a  log-normal  distribution  of  strength,  representing 
the  validity  in  ultimate  strength  of  the  material  itself. 

Therefore  extrapolation  has  been  done,  using  a  log¬ 
normal  distribution  with  mean  at  x  »  0  and  variance  se¬ 
lected  to  make  it  continuous  with  the  "t"  distribution  at 
x  =  -0.2 . 


This  leads  to  the  following  distribution  of  struc¬ 
tural  resistance : 


p(x)dx 


2a  a 

-IX- 


nlax2  +  x3)a 


-0.2  S  x  <  00 


.  (log10x+l)2 

4.6955  .002723 

(x+1) 


-1  £  X 


(3.6) 


s  -0.2 


with  a  =  0.05638 
x 

In  applying  the  reliability  function  it  is  often  an  advantage 
to  have  an  algebraic  expression  for  the  probability  distribu¬ 
tion  of  x  ('’’J  p(u).du^)  .  While  this  exists  for  the  "t"  dis¬ 


tribution  it  does  not  exist  for  the  log-normal  distribution 
and  a  polynomial  approximation  is  therefore  proposed  as  an 
alternative : 

10  3 

p(x)dx  =  0.6096  (x  +  1)  -1  £  x  *  -0.2  (3.7) 


-17- 


(c)  Probability  Distribution  for  Symmetric  Structures 


Some  aircraft  structures,  such  as  mainplanes,  are 
usually  symmetric  in  that  they  consist  of  two  nominally  iden¬ 
tical  halves.  Ultimate  failure  of  the  structure  then  consti¬ 
tutes  failure  of  the  weaker  of  two  members  from  the  popula¬ 
tion.  The  data  in  the  foregoing  analysis  are  from  asymmetric 
tests  or  from  symmetric  tests  on  asymmetric  structures  in 
which  failure  always  occurs  on  one  side,  which  is  also  nom¬ 
inally  the  weaker. 

Some  data  available  on  symmetric  structures  have 
been  analyzed  and  indicate  that  the  strength  of  the  two 
halves  is  not  independent  but  is  highly  correlated.  The 
preliminary  conclusion  is  that  it  is  more  accurate  to  apply 
the  foregoing  distribution  to  the  whole  structure  on  the 
assumption  of  complete  correlation  (i..e,. ,  both  halves  always 
of  identical  strength)  than  it  is  to  apply  the  distribution 
for  the  lesser  of  two  independent  ("t"  distributed)  variates 
(i.e.,  both  halves  selected  at  random  from  the  population). 

The  physical  significance  of  this  is  that  in  produc¬ 
tion  two  equally  high  (or  low)  strength  structures  are  pro¬ 
duced  at  the  same  time  and  assembled.  To  determine  the  pre¬ 
cis'-'  degree  of  correlation  (if  a  unique  value  exists)  will 
require  further  data  from  symmetric  structures. 
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The  effect  of  ultimate  failures  ir»  different  areas 


of  the  same  structure  is  represented  in  the  data  since  some 
types  of  structure  showed  both  tension  and  compression  fail¬ 
ures  at  the  ultimate  bending  moment.  However,  this  effect 
is  rather  small  in  uncracked  structures  since,  because  of 
general  yielding,  the  load  distribution  at  ultimate  failure 
is  characteristic  and  is  consistently  reproduced. 
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4.  ULTIMATE  STRENGTH  ANALYSIS 


Consider  a  load  S  chosen  at  random  from  the  popu¬ 
lation  of  service  loads  and  applied  to  a  structure  of  re¬ 
sistance  R  chosen  at  random  from  the  population  of  struc¬ 
tures.  Following  the  theory  developed  by  Freudenthal1  the 
probability  of  failure  is  given  by 

PF  =  II  Pr(R)  PS(S)  dRdS  (4.1) 

iRSS  } 

where  R  and  S  are  assumed  to  be  independent  and  the  inte¬ 
gration  is  carried  out  over  the  region  R  £  S  (,i.e . ,  the 
quadrant  bounded  by  R  =  0  and  R  -  S  as  shown  in  Fig.  1). 

(a)  Reliability  Equations 

Integrating  (4.1)  with  respect  to  R  and  S  over  the 
region  R  £  S 


CO 

CO 

- 1 

O 

4Pr 

(R)  ps(S)  dSdR 

00 

=  [ 

Pr(R) 

[1  -  Pg (R) ]  dR 

vo 

00 

■  I 

Pr(R) 

FS(R)  dR 

o 
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The  marginal  density  function  of  the  joint  distribution  (4.2) 
within  the  limits  0  <  R  <  S  is : 

p(R)  =  Pr<r)Fs(R) 

Multiplication  by  dR  produces  the  probability  that  failure 
will  occur  in  structures  with  resistance  in  the  interval 
R  to  R  +  dR.  This  is  evident  from  first  principles  since 
p  (R)dR  is  the  probability  of  structures  having  resistance 

Xv 

R  to  R  +  dR  and  F  (R)  is  the  probability  of  exceeding  a 

b 

service  load  R.  Hence,  assuming  independence 

pr  <R)  .dR  .F  (R)  = 

a  b  r 

where  P  denotes  failure  of  structures  with  resistance 
r 

R  to  R  +  dR, 

The  probability  of  failure  of  a  structure  with  re- 

A 

sistance  less  than  R  is  therefore 

A* 

R 

Pp(R<R)  =  f  pR(R)Fs(R) .dR  (4.3) 

j  o 

Changing  the  order  of  integration  Eq.  (4.1)  may  be  written 

to  R=S  a> 

Pp  =  J  as  J  ps(S).pR(R).dR  =  J  PS(S)  PR(S).dS  (4.4) 

"  o  J  o  "  o 

The  marginal  density  function  of  the  joint  distribution  (4.2) 
within  the  limits  R  <  S  <  ® 

P(S)  =  PS(S)  PR(S) 
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p(S).dS  =  Pr 

where  Pr  is  the  failure  of  structures  under  a  service  loads. 
This  also  depends  on  the  assumed  independence  of  R  and  S 
since  pR(S)dS  is  the  probability  of  occurrence  of  a  service 
load  between  S  and  S  +  dS  and  P  (S)  is  the  probability  of 

JR 

a  structure  with  resistance  less  than  S.  The  combination  of 

these  two  events  constitutes  failure  under  a  service  load 

between  S  and  S  +  dS .  Hence  the  failure  P  under  a  serv- 

r 

ice  load  S  may  be  expressed  as 

Pr  *  ps(S)  Pr(S) 

The  probability  of  failure  under  a  service  load  less 

/v 

than  S  is 

/N 

-  rS 

PF(S  <  S)  =  j  PS(S)  PR(S).dS  (4.5) 

J  o 

A  A 

The  relation  between  P  (R  <  R)  and  P  (S  <  S)  is  readily 
obtained  by  referring  to  Fig.  8  where  the  respective  areas 
of  integration  are  shown. 

Pp(R  <  R)  *  PF(S  <  h  +  Pr(R).Fs(R)  (4.6) 

which  is  obtained  by  the  following  consideration  of  first 


principles . 


A 

P  {failure  of  structure  with  R  s.  R}  ** 
r  1  J 

=  {failure  with  R  *  R  and  S  <.  R}  + 

P  {failure  with  R  £  R  and  S  <  R} 

r  L  J 

=  Pp  {S  <  R}  +  Pr  fR  SR}  •  Pr{S  >  R} 

=  Pp  {S  <  R}  +  PR  (R)  *Fg (R) 

Either  of  the  alternative  forms  for  P  in  Eqs. 

F 

(4.2)  and  (4.4)  may  be  adopted  depending  on  which  one  may 

be  more  easily  integrated.  The  functions  P  (R  <  R)  and 

F 

A 

P  (S  <  S)  are  of  interest  in  investigating  the  values  of 
F 

R  and  S  at  failure . 

Making  use  of  the  distribution  function  P  for 

3C 

the  relative  variation  in  structural  resistance  x  devel 


oped  in  Section  3  and  the  probability  distributions  F^,  F^ 
for  exceedance  of  service  loads  in  various  operations  listed 

A  A 

in  Section  2,  the  functions  P  {R  <  R}  and  P  {S  <  S]  may  be 

F  F 


J_1Pxtl£)fnt4nR)  •<* 
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■t 


U  or  an  may  be  expressed  in  terms  of  x  using  the  design 
R  R 

conditions  as  shown  in  Eq .  (4.2)  and  the  integral  may  then  be 


evaluated . 


A  5 

Considering  P  {S  <  S}  and  defining  y  =  * —  -  1 


with  Y  =  —  -  1  it  follows  that 


Pp{S  <  S}  =  Pp{y  <  Y]  =  J  ps(S)  PR(S).dS 


■  J  Vy)  t  px  (r- 11 -ds 


X 

=  f  p  (y)  Px(y)-dy 


(4.8) 


Y  dU 

=  if)  px(y)-dy 

Y  dn 

'  I_1fpnlnS)  Px(y,'dy 


dns  dus 

Here  again  ng  and  or  Ug  and  may  be  expressed  in 

terms  of  y  using  the  design  conditions  and  the  integral 
can  therefore  be  evaluated. 


!3ion 


In  formulating  the  design  conditions  it  is  assumed 
that:  (1)  The  distribution  of  the  applied  load  remains  con¬ 
stant  and  its  magnitude  is  proportional  to  the  acceleration 
produced  at  the  center  of  gravity;  and  (2)  Changes  in  air 
loads  due  to  gusts  are  proportional  to  the  gust  velocity. 
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The  following  notation  will  be  employed  in  the  dis¬ 


cussion  of  design  conditions: 


W  actual  operating  weight  under  different  conditions 

n  acceleration  at  the  center  of  gravity  due  to  the  ap¬ 
plied  load 

n  acceleration  at  the  center  of  gravity  in  straight  and 
level  flight  (normally  nQ  =  +1  for  positive  loading 
case  and  nQ  =  -1  for  negative  loading  case). 

An  acceleration  produced  by  a  change  in  the  applied  load, 

i  .e  . ,  n  =  n  +  An . 

-  o 

Kn  ^  (where  U  is  the  gust  velocity),  a  gust  sensitivity 
factor 


p  probability  that  a  specified  minimum  strength  R  is 
not  exceeded  in  the  population  of  structural  resist¬ 
ance  R,  i.e.,  p  =  P  {R  <  R  ). 

q  probability  that  a  specified  maximum  load  Sg  is  ex¬ 
ceeded  in  the  population  of  service  loads  S,  i.e., 

q  =  P  {S  ;>  S  } 
r  1  q J 

Rp 

aR  ,  the  ratio  of  the  specified  minimum  strength  to 

,p  ^R  the  mean  structural  resistance  ^ 

-  RP  qrVr 

v  a  safety  factor  defined  by  v  =  ~ 

P  #Q  b  O 

q  q 


S  has  been  taken  in  all  cases  as  the  value  speci- 

q 

fied  in  current  design  requirements.  Two  values  of  R^  have 
been  considered  in  order  to  apply  the  design  philosophy  of 


Reference  1 . 


mm 

To  simplify  the  notation  no  subscripts  are  used  on  v 

in  the  following  equations,  but  v  is  always  shown  associated 

with  a  corresponding  value  of  p  (or  a  ) . 

R  ,p 
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The  values  of  p  considered  are: 


(i)  p  =  0.1  based  on  the  assumption  that  Rp  is 

the  structural  resistance  corresponding  to  the  specified 

minimum  values  of  the  material  properties .  It  has  been 

9 

shown  by  Freudenthal  from  extensive  test  data  on  small 
specimens  that  approximately  10%  of  the  specimens  in  a  large 
sample  give  values  below  the  specification  minimum.  This 
indicates  that  a  value  of  p  =  0.1  is  reasonable;  the  cor¬ 
responding  value  of  a  based  on  the  "t"  distribution  in 

R,p 

(3.6)  is  a_  =  0.947. 

R  j  p 

(ii)  p  =  0.5  on  the  assumption  that  Rp  is  equal  to 

the  mean  structural  resistance  t±,(a  =  1.0).  This  is  as- 

R  R,p 

sumed  to  be  the  condition  relevant  to  present  aircraft  de¬ 
sign  standards,  since  in  current  design  the  ultimate  fail¬ 
ing  load  as  determined  in  the  usual  single  test  to  destruc¬ 
tion  exceeds  the  design  ultimate  load  by  a  negligible  margin. 

The  current  safety  factor  on  aircraft  design  loads 
is  1.5  and  hence  v=1.5,  p=0.3  is  the  equivalent  con¬ 
dition  for  comparison  with  service  conditions. 
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Considering  an  upward  load  applied  to  the  structure 

R  =  W  •  \  =  "("o  + 


Rp  =  *VVp"  vw  Cno  +  inq) 


R.a 


X  =  —  -  1  = 


R.P  .  x  =  (%  *  AnR)aR.P  _ 

R  (n  +  An  )  v 

p  o  q 


with 


AnR 


aR  .n 

2  0  and  hence  x  s  -Q  _  ]_  #  Hence 


(x  +  1)  (n  +  An  )  v 
_ 2 _ 2 _ 


(4.9) 


a. 


R,p 


n 

q 

a  .n 

-  n  ;  x  s  'g~ 9  -  1  (4.10) 

o  nqV 


or 


Anp 

UR  =  "k 

n 


(x  +  l)(n  +  K  U  ) v 
_ 2.  _  JLSLl 

a„  K 
R,p  n 


n 


_o 

C 

n 


„  ;  x  ^  a  .n 

K  R,p  o  _ 

*_0  '  1 


n 1  (4.11) 

Substitution  of  these  values  in  (4.7)  transforms 


F^AiIr)  or  Fu(Ur)  to  a  function  of  x  with  the  lower 
11  aR .  pn0 

limit  of  integration  taken  at  x  -  — Hr"9  -  1  or 

aR  p*no  °qV 

x  =  v 1 9  _  -  1.  For  a  downward  load  changing  the  sign  of 

K  U  v 
n  q  v 

n  produces 
o 


AHr 


(x  +  1)  ( An  -  n  )  v 
_ d  -  Q- 


a 


+  n 


(4.12) 


R,p 


or 


UR  * 


(x  +  1 )  (K  U  -  n  )? 
_ 0-2 - 2 _ 


n 


a  K 
R,p  n 


K 


(4.13) 


n 


with  aHr  and  UR  now  necessarily  greater  than  or  equal  to 
zero  for  x  i  -1 . 
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For  a  lateral  load  n  =0  and  therefore 

o 


AnR 


(x  +  l)Anq~ 


a. 


(4.14) 


R,p 


(x  +  1)U  ~ 


LL 


a. 


(4.15) 


R..p 


With  0  and  U  >  0  and  provided  that  x  >  -1, 

K  R 

<A 

Pp (R  <  R}  can  therefore  be  evaluated  from  (4.7)  by  sub¬ 
stituting  for  An  or  U  according  to  the  case  considered 

K  K 

Similarly  for  R  [S  <  S}  =  P  fy  <  Y]  with 


S  =  w(n  +  An  )  =  W(n  +  K  U  ) 
o  b  o  n  S 


(n  +  An  )a 

y  =  -  i  =  -2 - _  x 

^  (no  +  »q>v 


(4.16) 


by  comparison  of  (4.16)  and  (4.9)  it  follows  that  An  and 
AnS 

U  =  ~ —  are  obtained  directly  from  the  corresponding  ex- 
S  is. 

n 

pressions  for  An  and  U  by  replacing  x  by  y. 

R  R 

For  an  upward  load 

(y  +  1) (a  +  An  ) v  cl  n 

1  o  q  R.P  o 

An  =  -  - - -  n  ;  y  >  _ —  -  1 

S  aR,p  °  nqv  (4.0) 


(x  +  1) (n  +  K  U  )v  n 
U  = - 2 - n_2_  _  _o  . 

S  a. 


R,p 


K  ’  y  K  U  v 
n  q 


a  n 
R-lP-O  _  i 


(4.18) 
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1 


For  a  downward  load: 


Ans  = 


(y  +1) ( An  -  n  ) 
_ g  _  Q- 


a, 


+  n 


R,p 


(4.19) 


U, 


(y  +  1)(KU  -  n  ) v  n 

_  Q  g  -Q-  + 

K 


a_  K 

R,P  n 


n 


(4.20) 


For  a  lateral  load: 


An, 


(y  +  l)Anqy 


a, 


RiP 


(4.21) 


U, 


(y  +  l)Uqy 


a. 


R*P 


(4.22) 


(c)  Influence  of  the  Design  Parameters 

On  the  basis  of  the  preceding  equations  the  influence 
of  the  relevant  parameters  on  the  probability  of  failure  in 
civil  and  military  aircraft  structures  may  be  investigated, 
(i)  Civil  Transport  Aircraft 


Considering  Civil  Airworthiness  Requirements 


10 


the  feasible  range  of  the  parameters  may  be  established  as 


follows . 


For  positive  loads : 


Mw  *  (1  +  1.5) v  W 

K 


Maneuver  case 


“n  2  11  +Kn V5W  Gust  case 


For  negative  loads : 


M-r  *  (2  -  l)vW 


Maneuver  case 


2;  (K^U  -  1)^  W  Gust  case 
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Since  the  gust  loads  in  civil  aircraft  are  always 
relatively  severe,  the  gust  critical  condition  will  be  con¬ 
sidered  next.  Denoting  the  applied  load  in  the  steady 
flight  condition  by  n^W 


U_  =  a  R  >  (K  U  +  n  )  “  W 
R  R,p  p  n  q  o 

with  K  U  >1.5  for  upward  loads  (n  =  +  1) 
n  q  o 

K  U  >2  for  downward  loads  (n  =  -  1) 
n  q  o 

Minimum  values  are  thereby  established  for  the  de¬ 
sign  parameter  Kn .  The  influence  of  Kn  on  the  probability 

of  failure  is  indicated  by  U  ,  the  gust  velocity  corre- 

o 

sponding  to  the  average  structural  resistance  u  . 

R 

Setting  x  =  0  in  (4.11) 


vU 


U  =  - 

o  a 


R,p 


n 

K  ^a_  J 

n  v  R,p 


With  n  =  -  1,  U  increases  with  increasing  K  and  the 
o  o  n 


minimum  value  of  K  (K  =  77“  )  is  therefore  critical  for 

n  n  U 


q 

downward  loads  (since  v  >  1) .  with  n  =  +1,  U  decreases 

o  o 

with  increasing  and  the  upper  limit  of  is  critical 

for  failure  in  the  up-load  case. 

In  any  civil  aircraft  design  the  maximum  design  load 
is  never  likely  to  exceed  4W.  Considering  this  as  an  upper 
limit 
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(K^  +  l)vW  <4r<4vw 


•  •  K  U  <3 
n  q 


Realistic  values  for  the  design  parameter  K  can  therefore 

n 

be  expressed  as 

*  K  * 

U  n  U 

q  q 

Referring  to  (4.11),  the  effect  of  varying  the 

parameters  V  and  a  has  been  investigated  for  various 

R,p 

values  of  using  the  thunderstorm  gust  load  distribution 

from  Section  2a(i).  The  influence  of  K  has  been  investi- 

n 

gated  for  the  normal  value  of  v  -  1.5  with  p  =  0.1  and  0.5. 
For  comparison  a  calculation  for  a  particular  case  has  been 
made  using  the  gust  load  distribution  of  Section  2a(ii).  An 
outline  of  these  calculations  is  shown  in  Table  5  where  the 
values  of  the  parameters  are  listed  with  the  reliability 
equations . 


(ii)  Military  Aircraft 

The  probability  of  failure  under  gust  loads  is 
analogous  to  that  investigated  for  civil  transport  aircraft 
in  (i) .  However,  for  fighter  aircraft  the  maneuver  load 
spectrum  is  quite  distinct  from  the  exponential  gust  load 
spectrum. 
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The  semi-normal  probability  distribution  for  maneuver 
loads,  given  in  Section  2b(i),  is  not  readily  integrated  and 
the  form  of  the  reliability  equation  (4.8)  has  therefore  been 
adopted,  using  the  integrable  form  of  p(x)  in  (3.7)  with 
P( Ang)  . 

From  Eq.  (4.17) 

(y  +  1)  (n  +  An  )  v 
.  O  Q 

Anc  =  -  - 3 —  -  n 

S  a_  o 

R ,  P 

Hence  p(Ang)  may  be  transformed  into  a  function  of  y, 
introducing  the  parameters 


n  =  1  +  An '  ,  v,  a  and  a  with  n  =  +1 . 

q  q  R,p  n  o 


The  effect  of  the  specification  parameters  v  and  a  has 

R,p 

been  investigated  for  typical  values  of  n^  and  on*  The 

effect  of  the  service  load  parameter  an  has  also  been  in- 

2  1 

vestigated  taking  n^  =  8  —  and  7  —  corresponding  to  the 
FI  and  FII  categories  in  military  specification  MIL-A-8861 (ASG) 


for  the  standard  design  values  v=1.5,  a  =1.0.  The 

R,p 

cases  calculated  are  listed  with  the  relevant  values  of  the 


parameters  in  Table  4. 


(d)  Comparison  with  Service  Conditions 

To  test  the  application  of  the  foregoing  theory,  the 
probability  of  failure  has  been  predicted  for  specific  cases 


in  which  data  on  the  failure  rate  have  been  obtained  from 
service  records. 

(i)  CjLyil.i..TojQgBgr,l3 

Accident  statistics  for  United  States  and  United 

Kingdom  civil  transport  aircraft  in  the  post-war  period  have 

been  obtained  to  establish  an  overall  failure  rate  per 
11  12 

hour.  *  These  data  are  compared  in  Table  5  with  the 

results  predicted  by  theory  for  both  the  thunderstorm  gust 

data  and  the  general  gust  data  (Section  2a). 

Design  values  “=1.5  and  a  =1.0  have  been 

R,p 

taken  and  the  combined  failure  rates  under  up-gust  and  down- 
gust  are  presented  for  a  series  of  values  of  for  the 

thunderstorm  gust  spectrum. 

(ii)  Heavy  Bombers 

Data  from  a  large  fleet  of  heavy  bombers,  in  both 
high-  and  low-level  operations,  are  shown  in  Table  5.  Also 
presented  for  comparison  are  failure  rates  predicted  by  the 
reliability  equation  (3.7)  for  the  following  cases. 

(1)  Lateral  gusts  for  low-  and  high-level  opera¬ 
tion,  based  on  the  ultimate  gust  velocity  Uq  for  the  design 
operating  conditions 

(2)  Vertical  accelerations  (including  gusts  and 
maneuvers)  for  the  average  operating  conditions  of  the  fleet 
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at  both  high  and  low  level .  The  particular  service  load  data 
for  the  case  are  shown  in  Fig.  6. 

In  (1)  no  cases  of  wing  bending  failures  have 
been  reported  so  that  only  a  lower  limit  for  the  failure 
rate  can  be  established  from  the  service  operations. 

In  (2)  the  operating  conditions  at  both  high  and 

low  level  are  substantially  different  from  design  operating 

conditions.  The  value  of  U  taken  in  the  calculations, 

o 

although  based  on  experiment,  is  therefore  too  high  for  the 
low-level  operation,  and  too  low  for  the  high-level  condi¬ 
tions  . 

This  preliminary  estimate  of  the  failure  rate, 
therefore,  gives  no  more  than  an  approximate  indication. 

(iii)  Fighter  Aircraft 

Data  from  a  large  fleet  of  fighter  aircraft  in 
both  training  and  combat  operations  are  shown  in  Table  5.  In 
this  case  the  ultimate  load  at  failure  has  been  determined  on 
a  number  of  specimens  and  data  taken  on  V-g  recorders  during 
training  operations  are  also  available. 

The  failure  rate  has  been  predicted  by  the  relia¬ 
bility  equation  (4.8)  using  the  service  data  applying  to  the 
fleet . 
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5.  FATIGUE  STRENGTH  ANALYSIS 


Data  on  the  fatigue  performance  in  service  of  the 

three  classes  of  aircraft  being  considered  (.i.e^. ,  civil, 

bontoer  and  fighter)  are  included  in  Table  5.  These  have 

been  combined  with  whatever  fatigue  information  could  be 

obtained  on  the  structure  to  determine  the  risk  of  fatigue 

failure  r  (N) .  Following  the  proposals  in  Reference  1 

r  _.(N) 

tF  ~i 

f  =  - — |  has  been 

obtained  using  the  results  from  the  ultimate  strength 
analysis  in  Section  4. 

(a)  Civil  Aircraft 

Data  from  United  States  Air  Carrier  operations  over 
the  period  1946-1960  have  been  analyzed"^  and  are  presented 
in  Table  5.  Similar  data  for  British  aircraft  in  the  period 
1948-1956  have  been  reproduced  from  Reference  12.  The  total 
flying  hours  were  not  given  in  this  source  but  an  approxi¬ 
mate  figure  of  2  x  10^  hours  per  year  has  been  estimated 
indirectly  and  has  been  adopted  as  sufficiently  accurate 
for  the  present  purpose. 

The  average  failure  rate  per  hour  has  been  calcu¬ 
lated  for  both  fatigue  and  ultimate  load  failure.  The 
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United  States  and  British  figures  show  reasonable  agreement 
and  have  been  combined  to  give  an  average  estimate.  The 
calculation  of  a  uniform  failure  rate  is  valid  for  ultimate 
failure  since  it  is  regarded  as  due  to  random  causes.  The 
failure  rate  in  fatigue,  however,  increases  during  the  serv¬ 
ice  life  of  each  particular  aircraft  and  the  average  figure 
gives  only  a  gross  estimate.  It  has  been  multiplied  by  2 
on  the  assumption  of  a  linearly  increasing  rate  to  give  an 
approximate  figure.  At  the  very  low  probabilities  of  fail¬ 
ure  involved  the  risk  of  failure  is  virtually  the  same  as 
the  failure  rate. 

An  estimate  of  the  fatigue  sensitivity  of  (N*)  = 

PF  rF 

—  =  —  has  thus  been  made,  giving  the  overall  fatigue 
Pu  u 

sensitivity  for  the  whole  population  of  aircraft  at  the 
end  of  the  average  service  life. 

(b)  Heavy  Bomber  Aircraft 

A  sample  of  40  aircraft  in  the  fleet  has  been  con¬ 
sidered,  assuming  that  each  craft  was  operated  under  a  com¬ 
mon  spectrum  of  service  loads.  All  aircraft  were  inspected 
over  a  period  and  initial  fatigue  failures  were  noted.  At 
the  same  time  the  structure  was  modified  in  the  fatigue 
critical  area  whether  cracks  were  present  or  not. 
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From  the  data  listed  in  Table  6  it  can  be  seen  that 


19  aircraft  were  modified  before  initial  failure,  and  the 
data  therefore  include  19  arbitrarily  censored  values  of 
service  life.  Assuming  a  log-normal  distribution  of  life, 
the  mean  and  variance  have  been  estimated  from  the  follow¬ 
ing  maximum  likelihood  equations,  applying  to  a  sample  of 


n  with  n-k  censored  values: 


.15 


a  n-k 

hV  =  Z  +  ~  Z  A(v  ) 
Z  k  j-1  3 


(5.1) 


2  k  2 

oz  =  r  <z.  -  nz> 

i=l _ 


n-k 

k  -  Z  tv  A(v  )  ] 
j=l  3  3 


(5.2) 


where  Z.  =  log  N.  for  i  **  l,...,k. 

i  i 


2  .  —  (1 

v.  =  i  Z 


A(v.)  = 


-1/2  v2 
3 


f  e-1/2  dt 


v . 
3 


These  have  yielded  values  \i  -  3.2948,  q  -  0.005348,  giv- 

z  z 

ing  m-n  =  1970  hours  with  *  0.0731.  The  corresponding 
log  normal  distribution  has  been  drawn  in  Fig .  9  giving  the 
probability  of  failure  P(N)  at  N  hours.  The  risk 


-37- 


function  associated  with  P(N)  is 


rF(N>  = 


log  e  C e 


(log  N  -  kz) 

2az2  ) 


N-o 

z  *' 


e 

logN-p, 


-1/2  t  dt 
Z 


(5.3) 


A  plot  of  the  relative  frequency  of  the  individual  data 
points  has  also  been  made  by  employing  the  risk  function. 
Assume  that 


%• 


=  number  of  structures  remaining  at  life  N 

number  of  structures  in  the  original  population  which 
“  would  give  an  expected  number  m^  left  at  N 


Amn 


=  number  of  structures  that  fail  in  fatigue  during  a 

small  interval  N  to  N  +  iN  during  which  no  values 
are  censored  (M^  constant). 


The  model  postulated  is  that  censored  values  are  removed  at 


N,  and  no  censoring  takes  place  until  N  +  AN  when  further 


censored  values  are  suddenly  removed  and  so  on. 


L(N) 

mN 

=  5  « 

M„  is  not  known  from  the  data 
N 

p(N) 

ALUti 

an 

_  /  „ 

"  5T/ AN 

N 

r  (N) 

p(N) 

A\  /mN 

AmN 

L(N) 

Van/  “n 

m  .  AN 

N 

r(N)  is  therefore  the  only  function  independent  of 
Taking  the  relation  L(N)  =  e  io  e(t).dt 
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where 


where  N.  is  the  life  to  fatigue  failure  of  the  in  ,  the 
tr  tr 

structure  to  be  removed  from  the  population. 

The  frequency  distribution  for  the  21  values  of  fa¬ 
tigue  life  in  Table  6  has  been  plotted  in  Fig.  19  from  which 
it  will  be  seen  that  median  value  agrees  very  closely  with 
the  log  normal  distribution  but  that  the  variance  is  con¬ 
siderably  lower.  A  log  normal  distribution  with  a  standard  de¬ 
viation  s(logN)=.036  has  been  drawn  to  show  a  comparison 
with  the  frequency  distribution. 


-39- 


From  the  extensive  crack  propagation  data  avail¬ 
able  on  distributed  member  two-spar  structures,  a  crack 
propagation  period  of  60%  of  the  total  life  seems  realistic 
for  the  well-defined  initial  failure  considered  here.  This 
provides  an  approximate  estimate  of  5,000  hours  for  the  life 
to  final  failure. 

(c)  Fighter  Aircraft 

A  large  percentage  of  the  fighter  fleet  Was  inspected 
and  fatigue  cracks  of  varying  extent  were  found  in  a  number  of 
aircraft.  These  failure  data  have  been  corrected  on  the  basis 
of  known  crack  propagation  rates  to  correspond  to  a  typical 
initial  failure  for  all  cases. 

From  these  data  a  life  of  500  hours  with  a  standard 
deviation  of  0.186  is  obtained,  which  is  in  reasonable  agree¬ 
ment  with  test  results. 

There  is  no  record  of  complete  structural  failure  in 
service  due  to  fatigue  but  it  is  known  from  actual  test  data 
that  final  failure  occurs  in  another  part  of  the  structure. 

The  crack  propagation  time  is  approximately  20%  of  the  total 
life  and  it  is  therefore  not  surprising  that  no  initial  fail¬ 
ures  were  found  in  this  region.  The  fatigue  life  from  repre¬ 
sentative  test  data  is  3,100  hours  with  a  standard  deviation 


It  will  be  noted  that  the  standard  deviation  is  vir 

tually  the  same  as  that  for  the  bomber  and  is  in  agreement 

with  other  data  on  the  standard  deviation  at  relatively 

14 

short  fatigue  lives  as  quoted  by  Patching. 


(d)  Calculation  of  Fatigue  Sensitivity 


The  risk  function  r  (N) ,  given  in  (5.3)  can  be 

F 

transposed  into  a  non-dimensional  form  as  follows : 

.  o  £)2 


U' 


2a» 


V^N 


l0%/ 


e  .  e 


(N/V'°z  I 


e-1/2  tS.,t 


log  N/u 


N 


This  gives  aS  a  ^unct^on  —  which  applies  to 

N 

all  log-normal  distributions  of  the  same  variance  a  2 .  The 
function  has  been  plotted  for  a„  =  -073  in  Fig.  20  and 

Li 

used  to  calculate  the  fatigue  sensitivity  r  (N)  at  certain 
lives  as  shown  in  Table  6.  In  particular  the  life  Nq  for 
equal  risk  of  ultimate  and  fatigue  failure  has  been  estimated. 


6.  DISCUSSION  OF  RESULTS 


The  results  obtained  in  Sections  4  and  5  are  pre¬ 
sented  in  Tables  4  and  5  and  in  Figures  9  to  19. 

(a)  Ultimate  Failure 

In  Fig.  9  the  probability  of  failure  P  (X)  under 

F 

the  thunderstorm  gust  distribution  (Case  2  in  Table  4b)  is 
shown  for  various  forms  of  the  distribution  of  relative 
structural  resistance  x.  The  three  alternative  forms  for 
the  distribution  of  x  discussed  in  Section  3  (viz .  the  trun¬ 
cated  "t,"  the  truncated  "t"  with  polynomial  form  below 
x  =  -0.15,  and  the  truncated  "t"  with  log-normal  form  below 
x  =  -0.15)  show  good  agreement  for  the  total  probability 

of  failure  P„ .  However,  the  form  of  the  lower  tail  of  the 
F 

distribution  has  a  very  marked  effect  on  the  probability  of 

failure  as  is  shown  by  the  difference  between  the  high  order 

polynomial  and  log-normal  distributions. 

This  effect  is  also  pronounced  with  the  maneuver  load 

spectrum  as  shown  in  Fig.  16  where  the  probability  of  failure 

below  a  given  structural  resistance  P  (X)  and  the  probabil- 

£ 

ity  of  failure  below  a  given  load  P  (Y)  are  both  compared. 
In  fact,  for  most  cases  the  maneuver  load  spectrum  is  much 
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more  sensitive  in  this  respect  and  therefore  the  probabili¬ 
ties  of  failure  for  fighter  aircraft  have  been  truncated 
at  x  =  -0.25. 

The  curves  also  show  that  the  great  majority  of  struc¬ 
tural  failures  are  at  structural  resistances  near  the  mean. 
This  indicates  that  the  scatter  in  structural  resistance  at 
present  achieved  in  production  should  not  seriously  increase 
the  risk  of  failure  which  would,  therefore,  be  decisively 
determined  by  the  load  spectrum.  The  comparison  between  the 
total  probabilities  of  failure  predicted  by  the  invariate 
case  and  the  probability  distributions  substantiates  this 
view. 

However,  further  data  are  needed  to  investigate  the 
behavior  of  the  distribution  of  structural  resistance  at  the 
low  extreme  values.  It  should  also  be  noted  that  the  proba¬ 
bility  of  failure  is  very  sensitive  to  changes  in  the  mean. 
This  is  shown  by  the  curves  of  P  (X)  for  various  values  of 
the  design  parameters  v  and  p  in  Figs.  10  to  13  applying 
to  the  thunderstorm  gust  distribution. 

The  effect  is  clearly  shown  in  Fig.  15  where  the  total 
probability  of  failure  has  been  plotted  as  a  function  of  v 
for  the  thunderstorm  gust  and  maneuver  load  cases.  The  slope 
of  the  curves  near  v  =  1.5  shows  the  effect  of  a  variation 
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in  the  mean  structural  resistance.  Bearing  in  mind  that 
is  plotted  on  a  logarithmic  scale,  the  effect  is  very  marked 
particularly  for  the  maneuver  load  spectrum. 

This  suggests  that  the  probability  of  failure  in 
structures  developing  fatigue  cracks  may  be  very  much  greater 
than  anticipated.  However,  reliable  data  on  the  relationship 
between  structural  resistance  and  crack  length  in  the  general 
case  are  needed  to  investigate  this  important  question. 

The  influence  of  the  gust  sensitivity  is  shown 

in  Fig.  14  for  the  thunderstorm  gust  distribution.  The  ef¬ 
fect  on  the  probability  of  failure  is  not  aprreciable  except 
in  the  critical  downgust  case  at  low  values  of  K^.  However, 
there  is  a  lower  limit  of  Kn  =  .0303  applied  by  the  Air¬ 
worthiness  Requirements  in  this  case  as  discussed  in  Sec¬ 
tion  4  (b)  . 

For  fighter  aircraft  the  effect  of  the  mission  param¬ 
eter  a  is  shown  in  Fig.  17.  In  conjunction  with  the  val- 
n 

ues  of  an  listed  in  Table  1,  this  indicates  the  very  great 

influence  that  the  type  of  mission  has  on  the  probability  of 

failure.  In  changing  from  the  "general  mission"  role 

(a  =  1.27)  to  the  severe  ground  attack  (a  =  1.97)  the 
n  n 

probability  of  failure  per  load  application  is  increased  by 
five  orders  of  magnitude. 


In  Fig.  18  the  failure  rates  per  hour  of  the  three 
types  of  aircraft  (civil  transport,  bomber  and  fighter)  are 
compared  for  representative  missions. 

The  comparison  between  the  failure  rate  experienced 
in  service  with  that  predicted  for  average  operating  condi¬ 
tions  is  presented  in  Table  5(a). 

For  civil  transports,  an  average  cruising  speed  of 
200  mph  and  operating  height  of  20,000  feet  has  been  assumed 
representative  of  the  period  considered.  The  agreement  be¬ 
tween  the  failure  rates  per  hour  on  this  basis  is  very  good. 

The  comparison  for  the  bomber  in  the  case  of  verti¬ 
cal  acceleration  is  limited  by  the  relatively  low  total  hours 
accumulated  by  the  fleet.  The  failure  rate  predicted  for 
lateral  gusts  does  not  show  agreement  with  that  obtained 
from  service  records.  It  is  considered  that  this  is  mainly 
due  to  an  inadequate  estimate  at  this  stage  of  the  ultimate 
gust  velocity  Uq  corresponding  to  the  particular  operating 
conditions . 

The  failure  rate  predicted  for  the  fighter  is  based 
on  the  maximum  operating  conditions.  Since  the  operating 
weight  was  usually  less  than  the  maximum  allowable  a  failure 
rate  corresponding  to  r^  =  11  (all  upweight  «=  85%  maximum) 
has  also  been  presented.  The  comparison  of  the  latter  figure 
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with  the  failure  rate  recorded  in  service  is  quite  reason¬ 
able.  This  is  considered  a  realistic  comparison,  especially 
as  in  only  one  of  the  three  accidents  reported  was  the  take¬ 
off  weight  equal  to  that  allowable. 


(b)  Fatigue  Failure 

Fatigue  data  for  the  three  categories  —  civil  trans¬ 
port,  bomber  and  fighter  —  are  presented  in  Table  5(b). 

(i)  Civil  Transport 

No  data  referring  to  particular  aircraft  have 
been  analyzed  but  comparison  of  the  service  data  in  4(a) 
and  4(b)  indicates  that  the  ultimate  failure  rate  p 

u 

about  3/4  of  the  failure  rate  for  fatigue  p  at  the  end 

r 

* 

of  the  average  service  life  N  .  Considering  the  large  num¬ 
ber  of  hours  involved  and  the  close  agreement  between  the 
overall  figures  in  this  table,  the  factor  is  regarded  as 
significant . 

Accepting  the  assumptions  of  a  constant  value  of 

p  and  a  linearly  increasing  value  for  p  throughout  the 
u  F 

life,  this  infers  that 


PF(V 

PU 


rF(N  ) 

F  Q 

r  (N  ) 
u  o 


=  1 


at  3/4  of  the  average  service  life  N  . 
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(ii)  Heavy  Bomber 


The  fatigue  life  to  initial  failure  has  been 
assumed  to  have  a  log  normal  distribution  with  a  mean 

-  1,970  hours  and  a  standard  deviation  S(log  N)  =  .073 
N 

as  estimated  by  maximum  likelihood.  The  plot  of  the  ob¬ 
served  frequency  points  in  Fig.  19,  using  the  method  devel¬ 
oped  in  Section  5,  shows  good  agreement  with  the  median  of 
the  proposed  distribution.  However,  it  is  considered  that 
these  results  suffer  from  the  fact  that  the  inspection  in¬ 
tervals  have  become  effectively  shorter  as  the  total  hours 
of  the  fleet  increased.  Even  small  variations  in  the  extent 
of  the  detected  cracks  introduce  a  bias  under  this  procedure. 

The  rather  low  variance  for  an  initial  fatigue 
failure  obtained  by  maximum  likelihood  may  be  due  to  this 
effect  also,  but  the  overall  influence  on  the  mean  and  stand¬ 
ard  deviation  should  nevertheless  be  much  reduced. 

Calculation  of  the  fatigue  sensitivity  has  been 

carried  out  as  described  in  Section  5(d).  The  life  N 

o 

for  equal  risks  of  ultimate  and  initial  fatigue  failure  has 
been  estimated  as  960  hours.  This  value  is  relatively  short 
for  the  service  life  of  a  bomber  but  it  is  influenced  first 
by  the  fact  that  initial  fatigue  failure  has  been  consid¬ 
ered  and  second  by  the  very  severe  load  spectrum. 
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(iii)  Faaktgg. 


The  life  and  standard  deviation  given  in  Sec¬ 
tion  5  are  considered  to  be  very  well  based  and  since  the 
value  of  standard  deviation  is  appropriate  the  curve  in 
Fig.  20  has  again  been  used  to  deduce  the  values  of  fatigue 
sensitivity  shown  in  Table  6.  The  value  of  Nq  in  this 
case  is  1,900  hours  wihch  is  again  relatively  low,  espe¬ 
cially  considering  that  a  number  of  these  aircraft  had  ex¬ 
ceeded  a  service  life  of  2,000  hours. 
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7.  CONCLUSIONS 


(i)  The  reliability  theory  has  predicted  probabil¬ 
ities  of  ultinate  load  failure  in  good  agreement  with  those 
recorded  in  service  for  civil  transports,  a  typical  bomber, 
and  a  typical  fighter  aircraft. 

(ii)  The  variability  in  the  ultimate  strength  of 
structures  under  present  production  standards  is  satisfac¬ 
tory,  although  probably  marginal. 

(iii)  The  failure  rate  is  very  sensitive  to  the  mean 
structural  resistance  of  the  population,  particularly  under 
the  maneuver  load  spectrum. 

(iv)  In  view  of  (iii)  the  structural  reliability  of 
airframes  containing  even  small  fatigue  cracks  requires 
close  investigation  when  the  fail-safe  principle  is  invoked. 

(v)  The  failure  rate  is  greatly  influenced  by  the 
design  parameter  v  for  all  the  service  load  spectra  con¬ 
sidered.  Under  the  gust  load  spectrum  the  failure  rate 
increases  with  an  increase  in  the  gust  sensitivity  factor 
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* 

in  the  critical  downgust  case,  but  decreases  with  in¬ 
creasing  in  the  upgust  case. 

Under  the  maneuver  load  spectrum  a  variation  in  the 
failure  rate  of  5  orders  of  magnitude  is  predicted,  accord¬ 
ing  to  the  type  of  mission  considered. 

*  rF(N*> 

(vi)  The  fatigue  sensitivity  factor  f(N  )  =  ~ — 

u  v 

varies  greatly  according  to  the  value  taken  for  the  service 

life  N*  with  respect  to  the  mean  fatigue  life .  A  useful 

and  more  stable  indication  of  fatigue  sensitivity  is  the 

life  N  at  which  the  risks  of  ultimate  load  and  fatigue 
o 

failure  are  equal  [f(N  )  =  1]. 

o 

(vii)  For  airframes  in  the  three  categories  consid¬ 
ered,  civil  transport,  bomber  and  fighter,  the  risk  of  fa¬ 
tigue  failure  appreciably  exceeds  the  risk  of  ultimate  load 
failure  at  the  end  of  the  normal  service  life  N*. 
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Table  1 


Service  Load  Data 

FyfU)  »  Pr  [exceeding  gust  velocity  U(f.p.s)^. 

p  (An)  «  probability  density  function  for  acceleration 
n 

increment  An(^  ) 


r 

1.  Atmospheric  Gusts 

_ 

Type  of  Environment 

Total  gust 
counts  (+ve  and  -ve) 

(i) 

Thunderstorm  (Fig.  1) 

! 

FyOJ)  «  e‘°- 19713 

14  per  mile 

(ii) 

Civil  Airline  Routes  (Fig.  2) 

-0.344U  -0.208U 

FyOJ)  «  0 . 969e  ♦  0.031e 

i 

(see  Fig.  3) 

I 

,  s 

Ciii) 

Lateral  Gusts  (Fig.  4)  _Q  14JU 

Low  Level:  ^(U)  -  0.000345e 

j 

! 

I 

♦  0.«99655.-°-,68U  ♦  0.41e‘°'4H1U  ♦  0.09.-“°“ 

3 

2  x  10  per  hour 

» 

High  Level:  Fy(U)  -  0.032e"°'576U  +  0.968e"1,1SU 

i 

720  per  hour 

I 

2.  Maneuver  Accelerations 

Type  of  Operation 

Total  No.  of 
Acceleration  counts 

(i) 

Filter  Aircraft  (Fig.  5) 

?  .  &n)2 

Pnti")  m.$WFn  '  '  ^  "lth= 

1 

1 

1 

l 

i 

Ground  attack  A- A  (Fig.  5)0^“  1.97 

35  per  hour 

Ground  attack  C-C(Fig.  5)0^*  1.65 

40  per  hour 

Interceptor  F-F(Fig.  5)^“  1.70 

4<  per  hour 

General  mission  G“*  1.27 

n 

10  per  tour 

(11)  Heayy  Bowber  (with  equivalent  vertical  load 

factor  Anf)  (Fig.  6) 

An,  upward:  Ffcvn,)  ■  0.00314e  11,7  n*  ♦  0.997e  23,fo^nz 
An*  downward:  F(^'x)  ■  0.0022Se  ll*-39Anx 

N  per  hour 

♦  0.998e"28*78Anz 

N  per  hour 

with:  Anx»  (n^-0.72),  N  ■  6000  high  altitude  cruise 

Anx-  (0.68  -  nx),  N  -  6400  low  level 

Table  2 


Variation  in  Amount  of  Turbulence  with  Altitude 


Percent  flight  distance  in  turbulence 


Altit  ude , 


ft. 


Non-thunderstorm 
turbule nee 


Thunderstorm 
turbule nee 


0  to  10,000 
10,000  to  20,000 
20,000  to  30,000 
30,000  to  40,000 
40,000  to  50,000 
50,000  to  60,000 


18.0 

6.4 

4.5 
3.9 
3.4 
2.2 


0.10 

0.11 

0.062 

0.0067 

0.0017 

0 
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Tab  la  J. Data  on  Oltiaata  Strang th  of  Structuraa 
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Not  Significant  N.S.  10/20  10/20 
Significant  At  5*  2/20  3/20 
Significant  At  l£  8/20  7/20 


0363, .0-45-4 


-0.2  £  jr<- 


Calculate  P_  (y)  and  haaca  P_  (x)  for  tha  following 


C2*0--x)Ig-(«rt)JgS  2*6 
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i 


I 

1 

1  (iii)  Total 

c. 

i 

1 

i 

! 

1 

for  Mainplane  2  x  10 

0 

j<5  x  10‘7 

Combined 

. 

_ _ 

(i  and  ii)  2*6  x  10"9 

|  Filter 

1’06  x  105 

12-85  x  10~5 

Maneuver 

1 

1 

(Table  1  and  | 

Table  4)  ' 

n^-2,6^mi-97  3  x  10  ^ 

I  c+ 

r^-  11  I  6  x  10*  5  ! 


♦♦Refers  to  parti  _ulaf  type  of  bomber. 

^Calculated  for  cruising  speed  ■  200  raph ,  height  ■  20,000  ft. 

+ Based  on  measured  load  counts  in  service  *  S  counts  per  hour  with  a  >1.97. 

n 

(b)  Fatigue  Failure 


U.S.  Civil 

Air  Carrier 

U.K.  Civil 
Aircraft 

Total  for 

Civil 

Aircraft 

Hours 

No 

.of  Failures 

Failure 

rate 

(per  hour) 

Mean  Estimated  Life 

>0  x  106 

16  x  106 

66  x  106 

6 

7 

13 

2-4  x  10-7 

8-3  x  10”7 

4  x  10*7 

1 

1 

l 

i 

1 

j 

Heavy  Bomber 

r 

*  i 

(final  failure)!  2  x  10 

1 

10“  6 

5000  hours  (approx.) 

(initial 

A 

failure) 

2  x  10b 

19 

1970  hours 

Fighter 

(final  failure) 

1*06  x  105 

none de¬ 
tected 

3100  hours 

( initial 

i 

failure) 

1.06  x  103 

10 

480  hours  j 
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Table  6  -  Fatigue  Livet 
Initial  Failure  of  Heavy  Bomber 


N  1 

1  1 

Average  of  N 
for  interval  i 

Life  in  Hours  (N) 

- r 

— 

T  1 

,rjJ 

PCNp  «  1-e  Mti 

1  \i 

1298 

40 

* 

1434 

39 

1407 

1500  D 

38 

.0263 

.026 

i 

1504 

37 

! 

1529 

36 

i 

i 

i 

1538 

35 

1611 

34 

1708 

33 

1589 

1733  D 

32 

.0575 

.056 

1743 

31 

! 

1745 

30 

1743 

1753  D 

29 

.0919 

.088 

1754 

28 

1756 

1762  0 

27 

.1289 

.121 

1768 

1774  D 

26 

.1673 

.154 

1784 

1795  D 

25 

.2073 

.187 

'!  1812 

1830  D 

24 

.2489 

.220 

, 

1836 

23 

1 

1841 

22 

1340 

1855  D 

21 

.2965 

.257 

1856 

1858  D 

20 

.3465 

.293 

1364 

1870  D 

19 

.3991 

.329 

1877 

18 

! 

1877 

1884  D 

17 

.4579 

.367 

1898 

16 

1894 

1902  D 

15 

.5245 

.408 

1904 

1906  D 

14 

.5959 

.449 

1917 

13 

1924 

1949  D 

12  j 

.6792 

.493 

1994 

11 

i 

1984 

2010  D 

10 

.7792 

.541 

2019 

2028  0 

9  1 

.8903 

.590 

2030 

2031  D 

9 

1.0153 

.637 

2052 

7 

2052 

2075  D 

5 

1.1819 

.693 

2081 

2087  D 

5 

1.3819 

.749 

2107 

2128  0 

4 

1.6319 

.904 

2135 

3 

214C 

2158  D 

2 

2.1319 

.881 

2190 

1 

_ j 

D  =  Defect 
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Table  7  Estimation  of  Fatigue  Sensitivity  Factors 


(a)  Fatigue  sensitivity  factor  at  &even  life  times 


Life  to  Initial  Failure  i  Life  to  Final  Failure 


Y.-Q- 

• 

W  «  n  ; 

*■  , 

N  =  1/3  Mn 

• 

W  «  V^n 

_ | 

(i)  Bomber 

i 

/*n 

1970  hrs . 

1970  hrs.  j 

1 

SQ00  hrs. 

5000  hrs. 

ty  Mn  -9 

r.  ,nV  1.58  x  10 
(from  Fig.20)| 

-**  1 

7  x  10 

-9 

1.58  x  10 

-4 

7  x  10 

r. 

F 

-13 

8  x  10 

-7  i 

3.35  x  10 

-13 

3.16  x  10 

-7 

1.4  x  10 

r  (complete 
airframe) 

-7 

2  x  10 

! 

-7 

2  x  10 

-7 

2  x  10 

-7 

2  x  10 

from  Table  5 

i 

i 

! 

«  *V 

f(N  )  -  ~ 

ru 

-5 

4  x  10 

i.e 

-6 

1.6  x  10 

-1 

7  x  10 

N*= 

660  hrs. 

985  hrs. 

1670  hrs. 

2500  hrs. 

(ii)  Fighter 

i 

i 

/<n 

3100 

3100 

l 

r^n 

(from  Fig. 20) 

-9 

1.58  x  10 

-4 

7  x  10 

rf 

5.1  x  1013 

-7 

2.25  x  10 

r  (from  Table  5 
u 

) 

2.85  x  10'S 

2.85  x  10“5 

f(N*)  - 

■Lu 

1.8  x  10‘8 

7.9  x  10*3 

N*  - 

.1030  hrs. 

1550  hrs. 

I 

J 

-65- 


4PV 


PROBABILITY  OF  EXCEEDANCE  W  (  U  ) 


total  no  OF  GUSTS  (UP  a  DOWN)  per  nautical  mile 


Fig.  3.  Total  Gust  Counts  per  Mile.  General  Operation 
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Fig.  4.  Probability  Distribution  for  Lateral  Gusts. 


PROBABILITY  OF  EXCEEDANCE  F(nz) 


Fig.  6.  Vertical  Acceleration  Data  for  Heavy  Bomber. 
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Fig.  8.  Integration  of  Reliability  Function  in  Terms  of  R  and 


PROBABILITY  OF  FAILURE  PER  GUST  P- (X) 


10 


.-a 


10 


k-9 


10 


.-io 


10 


-1.0 


-0.6  -0.4  -0.2  0 

RELATIVE  STRUCTURAL  RESISTANCE  X 


Fig.  9.  Probability  of  Failure  in  Thunderstorm  Gusts. 
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PROBABILITY  OF  FAILURE  PER  GUST  P.(X) 


RELATIVE  STRUCTURAL  RESISTANCE  X 

Fig.  10.  Probability  of  Failure  under  a  Thunderstorm  Gust 
for  Various  Safety  Factors:  Upgust,  p  =  0.1. 
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PROBABILITY  OF  FAILURE  PER  GUST  PFIX) 


RELATIVE  STRUCTURAL  RESISTANCE  X 

Fig.  11.  Probability  of  Failure  under  a  Thunderstorm  Gust 
for  Various  Safety  Factors:  Upgust,  p  «  0.5. 
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PROBABILITY  OF  FAILURE  PER  GUST  PF  ( X ) 


RELATIVE  STRUCTURAL  RESISTANCE  X 

Fig.  12.  Probability  of  Failure  under  a  Thunderstorm  Gust 
for  Various  Safety  Factors:  Downgust,  p  =  0.1. 


TOTAL  PROBABILITY  of  FAILURE  PER  load  Pf 


Fig.  15.  Probabilities  of  Failure  as  a  Function  of  \. 
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PROBAB 


PROBABILITY  OF  FAILURE  PER  LOAD  Pf 


SERVICE  LOAD  PARAMETER  (7n 


Fig.  17.  Probability  of  Failure  of  Fighter  Aircraft  in 
Various  Missions. 
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PROBABILITY  OF  FAILURE  PER  HOUR 


RELATIVE  STRUCTURAL  RESISTANCE  X 


Fig  .  18. 


Comparison  of  Failure  Rates  under  Various  Load 
Spectra . 
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{n  >  h}  Jd  »  (N)d 
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Probability  Distribution  of  Fatigue  Life  for  Heavy  Bomber . 


NON  DIMENSIONAL  RISK  FUNCTION  fP(u*)fJ. 


